A generalized fluctuation relation for power-law distributions 
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Strong violations of existing fluctuation theorems may arise in nonequilibrium steady states char- 
acterized by distributions with power-law tails. The ratio of the probabilities of positive and nega- 
tive fluctuations of equal magnitude behaves in an anomalous nonmonotonic way [H. Touchette and 
E.G.D. Cohen, Phys. Rev. E 76, 020101(R) (2007)]. Here, we propose an alternative definition of 
fluctuation relation (FR) symmetry that, in the power-law regime, is characterized by a monotonic 
linear behavior. The proposal is consistent with a large deviation-like principle. As example, it is 
studied the fluctuations of the work done on a dragged particle immersed in a complex environ- 
ment able to induce power-law tails. When the environment is characterized by spatiotemporal 
temperature fluctuations, distributions arising in nonextensive statistical mechanics define the work 
statistics. In that situation, we find that the FR symmetry is solely defined by the average bath 
temperature. The case of a dragged particle subjected to a Levy noise is also analyzed in detail. 

PACS numbers: 05.70.Ln, 05.40.-a, 02.50.-r, 05.70.-a 



I. INTRODUCTION 

Fluctuation theorems have become a standard tool to 
characterize nonequilibrium states 0-0]. Independently 
of their specific formulation @~E1> a common underly- 
ing ingredient is an assertion about the symmetries of 
the fluctuations measures. For example, one of the most 
common formulation establishes that the probability of 
positive and negative fluctuations of a given variable dif- 
fer by an exponential weight proportional to the fluctua- 
tions magnitude. This FR symmetry has been confirmed 
in a wide class of experimental setups [l2| . While most of 
the analysis focus on variables such as entropy production 
and work (performed on a system), the FR symmetries 
may also be valid for nonthermodynamic variables [i~3| . 

A Brownian particle dragged by a spring through a 
thermal environment is one of the simpler arrangement 
where a FR symmetry can be theoretically predicted 
[3, [HI and measured [lfij]. While the work performed 
on the particle satisfy a standard or conventional FR, 
the heat fluctuations follow and extended FR 171. The 



exponential weight that relates the probability of positive 
and negative fluctuations does not scale linearly with the 
heat variable. Similar deviations with respect to a lin- 
ear dependence have been found in the injected power to 
systems driven by an external stochastic force [IH, [n| . 

Conventional or extended FRs [l7l - [l9j are essentially 
valid when the probability distributions obey a large de- 
viation principle (LDP) [2(| [2l|. This formalism also 
provides a solid basis for characterizing nonequilibrium 
states 22] . Therefore, stronger violations of conventional 
FRs should arise when power-law distributions [2l[ de- 
termine the fluctuations statistics. In fact, self-similar 
structures can not be studied in the context of a stan- 
dard large deviation theory [2(| |2l[ . 

One of the first analysis about the incompatibility of 
FRs with anomalous (power-law distributed) fluctuations 
was done by Touchette and Cohen in Ref. (23|. In 



contrast with previous Langevin models [lH . [l5j , where 
the environment influence is taken into account through 
a Gaussian white noise, they considered a stable Levy 
white noise. It was found that the ratio of the probabil- 
ities of positive and negative work fluctuations of equal 
magnitude behaves in an anomalous nonlinear way, devel- 
oping a convergence to one for large fluctuations. Hence, 
negative fluctuations of the work performed on the par- 
ticle are just as likely to happen as large positive work 
fluctuations of equal magnitude. This non- usual property 
strongly departs from the Gaussian case, where the va- 
lidity of a standard FR implies that positive fluctuations 
are exponentially more probable than negative ones. 

A similar analysis on anomalous fluctuation properties 
was performed by Chechkin and Klages [24| for the same 
kind of Langevin models. In the (power-law) Levy case 
the same conclusions were obtained. On the other hand, 
it was shown that standard FRs remain valid in presence 
of normal long time correlated fluctuations. The long- 
range correlations only alter the time-speed of the large 
deviation functions (LDFs) [201 ] ■ An analogous conclu- 
sion was obtained by Harris and Touchette [25| . 

In Ref. [26] Beck and Cohen introduced an alternative 
FR that arises by considering a superstatistical model 
[27l HH , where the particle environment develops (spa- 
ciotemporal) temperature fluctuations. As is well known 
[27I ] , this dynamics leads to power- law distributions aris- 
ing in Tsallis nonextensive statistical mechanics (29l - [35| . 
The alternative FR was derived by averaging tempera- 
ture in a standard FR. When the "measured" work dis- 
tribution develops power-law tails (26|, a very complex 
expression that does not has a clear physical meaning is 
obtained. 

From different points of view, all quoted analysis 
[2~il [26| sustain the main conclusion of Ref. 23] , that 
is, FRs in presence of power-law distributions acquire 
an anomalous (complicated) structure whose origin can 
be linked to the incompatibility of self-similar structures 
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with a (standard) LDP. The main goal of this paper is to 
introduce a generalized and alternative FR such that in 
presence of power-law tails the symmetry between posi- 
tive and negative fluctuations is expressed through a lin- 
ear dependence. Hence, a simple scheme for understand- 
ing previous results [23I4261 ] is obtained. Furthermore, we 
associate the proposed symmetry with a large deviation- 
like principle. In a long time regime, it allows expressing 
the FR through the symmetries of a set of LDFs associ- 
ated to the probability distribution and its characteristic 
function. As in the standard case [3| , both functions are 
related by a Legendre-Fenchel transformation (20j . 

Similarly to the case of standard FR, the proposed 
probability symmetry does not assume an underlying 
thermodynamic equilibrium, either extensive or nonex- 
tensive. Nevertheless, the alternative FR symmetry 
adopts a very simple structure when written in terms 
of a set of functions introduced in nonextensive entropy 
formalism [29j, |30j . The relevance of this stretched rela- 
tion becomes evident when considering previous related 
analysis (26|. Furthermore, it is expected that a FR as- 
sociated to a nonextensive (stochastic) thermodynamics 
may have a central role for understanding intrinsic fluc- 
tuations in nanoscale systems [361 ]. 

The paper is outlined as follows. In Sec. II we intro- 
duce the generalized FR symmetry. In Sec. Ill we apply 
the alternative definition to a specific physical system 
characterized by power-law tails. It is analyzed the fluc- 
tuations of the work done on a dragged particle immersed 
in a complex environment that develops (spatiotemporal) 
temperature fluctuations 1261 . Furthermore, the Levy 
model introduced in Ref. [23[ is analyzed in detail. In 
Sec. IV we give the conclusions. In Appendix A we 
develop a large deviation-like principle and characterize 
the symmetries of the LDFs. These last results rely on a 
saddle-point approximation performed in Appendix B. 
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where now £ p h play the role of a (physical) conjugate 
variable. Notice that even with this parameter redefi- 
nition, the proposed symmetry can be written solely in 
terms of p(x). The generalized FR written in terms of £ 
is called "unnormalized scheme," while in terms of C, p h 
"normalized scheme." 

The previous relations, Eqs. (JT|) and @, strongly de- 
parts from the standard ones. Their structure is sim- 
plified by writing them in terms of a set of functions 
introduced previously in the context of Tsallis nonex- 
tensive statistical mechanics [29]. A q- logarithm and q- 
exponential functions are defined respectively as 



\n q x 



A-q _ i 
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[1 + (1 - q)x] i-« . 



(4) 
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jointly with the generalized q-product operation 

x® q y = [x^ + y 1 ^- 1]^?. (5) 

The relevance of these definitions comes from the prop- 
erty that in the limit q — > 1 they recover the stan- 
dard logarithm [lni x = In a;], the exponential function 
[exp-L x = exp x] , and respectively the standard product 
[x ®i y = xy\. In terms of them, Eq. ([1} can be written 
as 

ln q \p(x)} - ln q [p(-x)} = (x, (6) 
while Eq. © becomes equivalent to 

p(-ar) = p(x)® q e q Cx , p(x) = p(-x) ® q e q x . (7) 

Hence, the proposed FR can be read as a "deformation" 
of the standard ones. The parameter q measures the 
degree of departure with respect to the standard case 
(q — > 1). Notice that by using the properties \n q (x<S> q y) 



e v [29(, the consistence 



*g ^ q 



For an arbitrary stochastic variable x s t, with proba- 
bility distribution p(x), a standard FR symmetry is de- 
fined by the relation ln[p(x) /p(—x)] = £x, or equivalently 
p(±x) = p(^fx)e^ x , where £ is a positive constant. Here, 
we propose the alternative FR symmetry 



\p(x)] 



1-9 



\p(-x)] 



1-9 



(x, 



i-q 

which in turn can be expressed as 

p(±x) = {\p(Tx)] 1 -«±(l-q)(x}Th 



(1) 



(2) 



between the previous two expressions becomes evident. 

The standard FR symmetry can also be writ- 
ten in terms of the characteristic function Z(X) = 

fZo dx P( x ) e ~ XX ' that is > Z ( A ) = Z (~ X + 0, whid i in 
turn implies Z(C) — 1. Taking into account Eq. ([2]), it 
follows that here we must consider the generalized ex- 
pression 

/>oo 

Z q (X)= / dxMx)} 1 -* - (1 - q)\x}^, (8) 



(9) 



which can be rewritten as 

Z q {\) = / dx[p(x) «), e q Xx ] 



The real parameter q is related with the index of the 
power-law tails. Its domain will be specified later on. 
The units of the constant £ are [l/a;] 2-9 . Nevertheless, in 
the standard case it is a conjugate variable of x, that is, 
its units are [l/a;]. This discrepancy can be avoided by 



By usin g th e associative property x < 



y) ®q 



|, the generalized FR [Eq. 



>q (V ®9 Z ) = ( X ®? 

© or 0] implies 



the equivalent symmetry 

z q (\) = z q (-x + o, 



z q (0 



(10) 
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Notice that this condition is exactly the same that defines 
the standard case, q — > 1. 

The expression ([9]) naturally arises in nonextensive sta- 
tistical mechanics. Hence, it properties are well known 
[29j . Z q (X) is the generating function of a kind of gener- 
alized moments of x s t , which are calculated from powers 
of p(x). 

In accordance with the normalized scheme [Eq. (|3"|)]. 
we introduce the definitions 



Zf{X ph ) = Z q (X(X p h)), A 



A 



ph 



(11) 



Thus, Eqs. Q and (|T0|) allows us to write the equivalent 
symmetry 

Zf(X ph ) = zf(-x ph +c P h), zf(c ph ) = 1. (12) 

A. q-Gaussian distributions 

The relations ([6]) and (|7|) define the generalized 
FR. Equivalently, it can be expressed through the "q- 
characteristic function" ([9]) leading to Eq. (fT0|) . Here, 
we search which kind of distributions may satisfy these 
relations in an exact way for any value of x. 

Gaussian distributions always satisfy the standard 
symmetry corresponding to lim q — > 1 . Normal distribu- 
tions emerge naturally when formulating a central limit 
theorem, as solutions of linear Fokker-Planck equations, 
or by maximizing Gibbs entropy under a second moment 
constraint. Similarly, q-Gaussian distributions 1291 are 
related to a generalized central limit theorem [31| . are 
solutions of a kind of non-linear Fokker-Planck equations 
[32| . and maximize nonextensive Tsallis entropy under a 
generalized second moment constraint [33j |. They read 



p{x) 



VP 



exp J-/3(x - x y 



(13) 



where /3 _1 measures the width of the distribution and 
Af q is a normalization factor such that J_°° p(x)dx = 1, 
(q < 3). For 1 < q < 3, p(x) is characterized by power- 
law tails, p(x) ~ l/(x — xq) 2 /^ 1 }, which is the case of 
interest in this paper. Hence, the index q determines 
the exponent of the power-law tails. The first moment 
of p(x) is finite for 1 < q < 2, while the second one for 
1 < q < 5/3. In the domain 1 < q < 3, the normalization 
constant reads fe9J 



,, _ ; 7T r (2( 3 g-l)) 

9 ~^-i rc-ij) : 



1< q < 3, 



(14) 



where T(y) is the Gamma function. 

It is immediate to prove that the generalized FR ^ [or 
Eq. |J7|)] is satisfied, for any value of x, by the q-Gaussian 
distribution (| 13[) with 

Q 3 ~ q O 

C = ^j2_, Cph= 4x p(-^-). (15) 



In the second expression we used the integral 



+ OG 



(16) 



valid for the distribution (fTB")) . When g — >■ 1, it follows 
the standard Gaussian expression £ = Q p h = Ax$f3. On 
the other hand, it is simple to realize that an arbitrary 
distribution p(x) satisfy the symmetry ([6]) when x is re- 
stricted to the interval where it develops power-law tails, 
p{x) ~ l/(x — rEo) 2 ^ 9-1 ' • Hence, even in presence of 
power-law tails, a linear relation characterize the sym- 
metry between positive and negative fluctuations of equal 
magnitude. 

Using that J_™ exp g (— ax 2 ) — Af q /^/a, after some al- 
gebra the q-characteristic function © associated to the 
q-Gaussian distribution (fl3|) reads 



z q {x) = exp- {(: 



A 2 AC"« 



4/3^ 



2 AytfJ 

where the index of the exponential reads 

1 + 5 



Xx 



q = 



3-q 



(17) 



(18) 



Consistently, the function (fT?)) satisfies the symmetry 
PU|) with the constant C given by Eq. (|T5|) . On the 
other hand, Z q lh {X) reads 



\ 2 



4/3(^2) 



X p hXo 



(19) 



B. Fluctuation relations for time-scaled variables 

While in the previous proposal we did not include time 
as an explicit parameter, its generalization to time depen- 
dent variables is immediate, x s t — > x st (t), p(x) — ¥ p(x, t). 
On the other hand, when studying nonequilibrium sys- 
tems it is common to define the variable of interest as a 
time-scaled one, 



1 f* 
Hst{t) x — / dtv s 
1 Jq 



(t) 



x st (t) 



t K 



(20) 



where (d/dt)x st (t) — v st (t). Usually, the time scaling is 
proportional to the elapsed time, k = 1. Hence, fi s t(t) 
can be read as a "time-average" velocity. Here, we adopt 
a more general point of view by considering arbitrary 
values of the exponent n > 0. In general, the definition 
(|20p makes sense in an asymptotic time regime. From 
now on, we use the symbol x for denoting an equality 
valid in a long time regime [20]. Clearly, the conditions 
that guarantees the achievement of this regime depend 
on each particular system. 

For variables such as fi s t(t) one can also define a gen- 
eralized FR. While its structure is very similar to the 
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previous case, it is worthwhile to write it explicitly. The 
probability distribution p(/J>) of p> s t{t) follows from the 
change of variable p(fi)dfi = p(x)dx. Taking into account 
that p(x) satisfies Eq. ©, a natural extension of the 
generalized FR is 



{\n q [p(n)} - ln 9 [p(-/i)]} x afi, 



(21) 



or equivalently 

p(-ri^)«,e^, P(j*)~p(-I*)® 9 e°i*\ (22) 

The constant a as well as the exponent rj > depend on 
each particular problem. Similarly the q-characteristic 
function of fj, s t{t) is defined as 

poo 

Z q (X) = / dn[p(n)® q e- xt ^}. (23) 



From Eq. (|2"2"|) , it satisfies the symmetry 



Z q {\) x Z q (-X + a), 



Z q (a) x 1. 



(24) 



Notice that in order to simplify the notation we do not 
explicit either that Z q corresponds to the variable /i s t{t) 
or its dependence on time. 

For the normalized scheme, Eq. ([3]), we define the 
equivalent FR 



i{H[p( M )]-lnX- M )]}- 



A*. (25) 



In this case, the normalization time factor is defined with 
a different exponent, S > 0. By introducing the charac- 
teristic function ZP h (X p h) = Z q (\(\ p h)), where 



A 



A 



ph 



it follows the equivalent symmetry 



(26) 



(27) 



When the probability distribution p(p) has an asymp- 
totic (t — > oo) exponential structure (q — > 1), the (stan- 
dard) FR can be analyzed through a large deviation 
theory [H HJ. The LDFs, that is, the factors that 
scales the time dependence of the probability distribu- 
tions and the characteristic functions, are related by a 
Legendre-Fenchel transform. Furthermore, the relation 
([2Tj) (<?—>• 1), or equivalently (|25|) . implies some symme- 
tries for the LDFs [J]. In the Appendixes, by establishing 
a generalized large deviation-like principle, we demon- 
strate that these results can be generalized to the present 
context (q ^ 1). Furthermore, general relations linking 
the unnormalized [Eq. (|2ip ] and normalized schemes [Eq. 
pj]) ] are formulated [Eqs. and (|A"27| ]. 



III. WORK PERFORMED ON A PARTICLE 
IMMERSED IN A COMPLEX ENVIRONMENT 

In the previous Section we have defined the alternative 
FR and in the Appendixes we developed a related large 
deviation-like principle. Here, we apply the proposal to 
a specific physical example. We study the fluctuations of 
the work done on a dragg ed particle [U, HH immersed 
in a complex environment able to induce power-law tails 
in the particle statistics 0, [26| . While this model has 
been studied previously, the following analysis show the 
conceptual and physical relevance of proposing an alter- 
native FR consistent with a LDP. 

The position x(t) and velocity v(t) of the particle obey 
the equations 



dx(t) 

dt 
dv(t) 

l : 

dt 



v(t), (28a) 
-jv(t) - k[x(t) ~ x*(t)} + £(t). (28b) 



The contribution —jv(t) gives the damping force. The 
term — k[x(t) — x*(t)} is the force produced by an har- 
monic potential where the position of its minimum is 
given by the arbitrary function x*(t). Finally, the envi- 
ronment influence is introduced through the noise 

In the overdamped regime, mk <C 7 2 , the position evo- 
lution can be approximated by the stochastic equation 

7^ = -*[*(*)-**(*)] +£(*)• (29) 

The work performed on t he p article by the harmonic drag 
force during a time r is |l5l. l23j 



[x{t)-x*{t)]v*{t)dt, 



(30) 



where v*(t) = (d/dt)x*(t). 

Independently of the environment model, the noise 
probability distribution is symmetric around the origin. 
Hence, when it exists, the average noise intensity is null, 
(£(*)) — 0- By assuming valid the same property for the 
initial particle position, (x(0)) = 0, we introduce the 
function 

M T = kf dtv*(t) f dt'e-^-^/^v*^'), (31) 
Jo Jo 



where the characteristic time is 

_ 7 
T °-k- 



(32) 



If the distribution of the noise admits to define a first 
moment, it is simple to realize that M T gives the average 
performed work, M T = (W T ) . In general, it only defines 
the most probable value of W T . 

For the Langevin dynamics (|29|) . transient behaviors 
develops when r < Tq. Therefore, the long time regime 
is achieved when t > ro. Under this condition, the ex- 
ponential factor in Eq. (1311) can be approximated by a 
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delta-Dirac function, e * x To<5(t — t'). Hence, we 
get 



M T x 7 / <%*(t)] 2 , 
Jo 



(33) 



where we have used that fcro = 7. 

The object of interest is the probability density P T (W) 
of performing a work W up to time r. We also study 
the asymptotic statistic of the dimensionless stochastic 
variable 



(34) 



whose probability density is denoted as p T (w). At any 
time, the most probable value of w T is one. Equivalently, 
when it exists, (w T ) = 1. Notice that W T and w T cor- 
respond respectively to the variables x st and /i st of the 
previous section. 



A. Standard model 

In order to clarify the next results, here we briefly re- 
view the standard case of a thermal environment at tem- 
perature T. Therefore, is a Gaussian noise whose cor- 
relation, consistently with a fluctuation-dissipation theo- 
rem, is (f (t)f (s)) = 2p- 1 1 8{t-s), where = k B T. As 
demonstrated in Ref. [15| . for thermalized initial condi- 
tions, P T (W) is a Gaussian distribution 



Pr(W) = 



1 



1 exp 



(W ~ M T ) 



2V T 



V^rVr 

where the time dependent width V T is 

V r = 2/3~ 1 M T . 
P T (W) satisfies the standard FR symmetry 
ln[P T (W)] - ]n[P T (-W)} = PW. 



The proportionality with /3 arises owe to the relation I 

After a simple change of variables, the probability den- 
sity of {3D is 



(35) 



(36) 



(37) 



p T {w) 



1 



y/2irV r 



exp 



(w - w ) 2 



2V T 



where 

W Q = 1, 

Hence, it follows the symmetry 



T M? (3M T ' 



^j{\n\p T (w)] - \n[ P r{-w)}} = ^-w 



(38) 



(39) 



(40) 



The exponent 8 — r\ is chosen in such a way that in 
the asymptotic regime the contribution M t /t s does not 



depends on time. We assume an accelerated potential 
movement 



x*{t) = 



V*t L 



(1 + a c ) ' 



v *(t) = vX c , 



(41) 



where u* is an appropriate constant and a c > 0. The 
usual assumption of constant velocity, v*(t) = u„, is re- 
covered with a c = 0. The average work (|3"3")l behaves as 



M T 



1 V * _ T l+2a c 



1 + 2a c 

Therefore, the FR gD]) is defined with 

1 



8 = 1 + 2a c , 



*ph 



1 + 2a, 



(42) 



(43) 



Notice that in general 8 7^ 1. A similar situation emerges 
in p resence of normal long time correlated fluctuations 
[24l [25j. Here, this dependence arises from the accelera- 
tion of the harmonic potential, Eq. (|4"Tj) . In fact, 8=1 
when the velocity is constant, a c = 0. On the other hand, 
owe to the chosen normalization, Eq. (|34p . a p h is not only 
proportional to /3 (factor 7U*). 



B. Superstatistical model 

Superstatistics [UHil consists of superpositions of dif- 
ferent statistics for driven nonequilibrium system with 
spatiotemporal inhomogeneities of an intensive param- 
eter, such as for example the inverse temperature /3 in 
the previous example. We assume that the time scale on 
which /3 fluctuates is much larger than the typical fluctu- 
ations of W T , that is, To = j/k. Hence, the work distri- 
bution can be written as P T {W) x / d/3/(/3)P^(W, (3), 
where P^ (W, f3) is the probability density of the normal 
environment case, Eq. (|35|) . This limit was analyzed in 
Ref. [26| without resorting on a LDP. A very complex 
probability relation arises, having not any clear depen- 
dence with the average environment temperature. These 
drawbacks are surpassed with the present approach. 

As in Ref. [27J, the temperature fluctuations are de- 
scribed by a Gamma distribution 



fifi) 



1 



■(— 

r(n/2) V2/3 



i/2 



a-1 



exp 



( _ 

V 2/3 ^' 



(44) 



where /3o corresponds to the average inverse temperature, 
Po = J d(3f((3)(3, and n is a positive constant. After 
averaging over this distribution, it follows the q-Gaussian 
distribution 



Pr(W) 



1 



2N*V T 



: exp 



(W - M T ) 2 
2V T 



Here, the width function V T is given by 
V T = 2/3" 1 M r . 



(45) 



(46) 
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Notice the similitude with Eq. (I36|) . Here, M T also fol- 
lows from Eq. ([551) . On the other hand, the characteristic 
parameters are [22| 



1 



1' 



3-g 



Ad- 



(47) 



It is difficult to extract some physical information by 
analyzing the distribution (|45|) through a standard FR 
symmetry [2f3], Eq. ((3Z|) In fact, given that P T (W) is 
a q-Gaussian distribution, it satisfies the generalized FR 
©. The parameter £, from Eqs. (fTS"]) and (|4"5")) can be 
written as 



c 



q-\ 3-g 



Cph = /?0 



(48) 



While the unnormalized FR is defined by a linear depen- 
dence on W, C depends on time. Hence, its physical con- 
tent is unclear. Nevertheless, in the normalized scheme, 
Qph corresponds to the physical (average) temperature 
/3o- Explicitly, the work probability distribution satisfies 
the FR symmetry 

Z T q {\n q [P T (W)} - \n q [P T {-W)]} x /W, (49) 

where for notational convenience we have defined 



7 T 



dW[P T {W)] q . 



(50) 



As in the standard case, in the present approach the envi- 
ronment (average) temperature is the scaling parameter 
of the probabilities linear relation (14"9"1) . Furthermore, this 
novel and elegant result recovers the standard FR (|37|) 
in the limit q — > 1, that is, in absence of temperature 
fluctuations, f{fi) = S(/3 - /3 ). 

The probability density of w T , Eq. is 



p T (w) 



1 



2M*V T 



: exp„ 



(w - w f 



2V T 



where the coefficients read 



w 



1, 



(51) 



(52) 



M T 2 (3 q M T ' 
Hence, p T {w) satisfies the symmetry 

-^{ln q \p T (w)\ - \n q [p T (-w)}} x w X a ph w, (53) 

where the normalization reads 

dw[p T (w)] q . (54) 

This integral can be performed by using Eq. fj 16|) . By 
assuming the velocity dependence (|4ip . it follows 



1 + 2a c , 



1 



a P h 



1 + 2a, 



(55) 
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FIG. 1: (Color Online) (a) Probability distribution p T (w) for 
the superstatistical model, Eq. (|51|) . for different times r. 
The inset details the peak region, (b) Plot of Zg {ln q [p T (w)] — 
ln q [p T (—w)]} for each time r. The inset shows the collapse of 
these curves to a single line when plotting the normalized FR 
symmetry $(w) = zl {\n q [p T (w)]-ln q [p T (-w)]}/T S , Eq. (1551) . 



The potential velocity is constant, a c = 0, Eq. (|41[) . Hence, 
5 = 1, Eq. (|55|l . In the natural units of the problem (see text) 
the parameters are v* = 7 = k = 1. The average temperature 
is Poivl/k = a ph To = 0.3, and q = 1.8, [Ea.(|i7|]. 



The previous analysis demonstrate that all features of 
a standard FR remain valid for the present example if one 
use the generalized FR. This result is explicitly shown in 
Fig. 1. 

In Fig. 1(a) we plot the q-Gaussian distribution (jSTj) 
for different times r. In the inset we show the peaks of 
the distributions. Their wide diminish with time. We 
assumed a constant potential's velocity; a c = in Eq. 
(|4"Tj) . Hence, 6 = 1. By using the (natural) units of mass 
too = j 2 /k, distance xo — v*j/k, and time tq = 7/fc, it 
follows 7 = k = u* = 1. Therefore, the unique free pa- 
rameters are q [Eq. (|4?p] and the (dimensionless) noise 
intensity, that is, the average temperature of the distri- 
bution (gU), /3 (7«*) 2 /fc = a ph To. 

In Fig. 1(b) we plot the dependence with w of 
z^{\a q [p T (w)] — ln q [p T (— w)]}. The index g is the same 
that defines the q-Gaussian distribution of Fig. 1(a). For 
each time, a linear behavior is evident. In the inset we 
show the collapse to a single line when introducing the 
time normalization factor (1/r 5 ), Eq. ()53|) . 
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In the unnormalized scheme [Eq. (|2"T1) ]. the FR reads 



— {\n q [p T {w)] - ]n q \p T (-w)]} 



where the coefficients are 

,,2 . Ina 



\P q M T ] 



T-n{2N q ) 1 -i 



aw, 



1 



a 



(2Af q y-i \l + 2a. 



3-9 



(56) 



(57) 



Here, the coefhcient a does not have a clear physical 
meaning. Nevertheless, as commented before, based on 
a large deviation-like theory (Appendixes) it is possible 
to establish some general relations between a and a p h, 
as well as between r\ and 8 [Eqs. (|A3|) and (|A27|) ]. which 
are satisfied in the present case. 



C. Levy noise model 

In Ref. [23[ the noise £(£) was taken as a symmetric 
stable Levy noise. Furthermore, different experimental 
setups where the model may be explicitly measured were 
proposed. While the analysis presented in that contribu- 
tion is completely right, here we study the same problem 
by using the generalized FR. We explicitly show that a 
FR can be established only when the probabilities satisfy 
a LDP. 

The noise is defined by its characteristic functional 

G e [X(t)] = (expi J dt\(t)£(t)y (58a) 

= exp(-6^ cft|A(t)| CT ), (58b) 

where X(t) is an arbitrary test function. The constant b 
measures the noise intensity and < a < 2. Due to the 
linearity of the stochastic dynamics (|29|) . the work (|30|) 
is also a stable variable with the same index a. It char- 
acteristic function, G^(A) 
then reads 



J+™dWP T (W)eiq>[i\W], 



G W (X) = (exp[iAW T ]), 

= exp(iAf T A - B T |A| C 



(59a) 
(59b) 



where M T is defined by Eq. ([31]) and B T can be obtained 
after writing Gw{X) in terms of Gf [A(i)]. For arbitrary 
velocities v*(t), we get 



B 



k\< T 



dt 



dt'e-^-W^v^t') 



(60) 



For r ^ To, M T can be calculated from Eq. ([3"3")l . while 
B T: after taking e - '* _t )/ T ° x T S(t' — t), can be approx- 
imated as 



B T -b f dt\v*(t)\ a 
Jo 



(61) 



Eq. (|59|) corresponds to the Fourier transform of a 
Levy probability distribution. As is well known [37| . for 

\W -M T \> B\ ja it develops power-law tails, 



Pt(W) 



c a B T 



|W-M T | 1 + CT ' 



(62) 



where c a = ir 1 a sin(7rcr/2)r(cr). Only when a = I, one 
gets a simple analytical expression valid for any value of 
W, P T {W) = (V T /n)[(W - M T f + V?}- 1 . It is expected 
that Eq. ([62]) satisfies the normalized FR 



Z r q {\n q [P T (W)} - ln g [P T (-W)]} ex C P h(r)W, (63) 

where the symbol ex denotes both an asymptotic time 
regime (r > r ) and \W - M T \ > B\ la that is, values 
of W in the power-law regime. The parameter q and the 
function £ p h (t) can be found by mapping the approxima- 
tion (1621) with the power-law behavior of the q-Gaussian 
distribution (|45|) . We get 



V T = C'B< 



(64) 



where c' a = ^(j^) 1 ^ (c a J\f q )i . Notice that 5/3 < q < 3 
jiT'li . With these relations at hand, the time dependent 
function C(r) reads 
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with c" 



(65) 



In deriving this expression we 



c' a l+a l+2a c 

assumed the general velocity dependence (|4"Tj) . While in 
the power-law regime the Levy distribution satisfy the 
generalized FR (|6"3"j) . the proportionality constant [Eq. 
(|4*5|] becomes time dependent, ( p h — > Cp^( t )- At long 
times, for any value of a 6 (0,2), it vanishes. Hence, 
consistently with the results of Ref. [23], one conclude 
that asymptotically positive and negative fluctuations of 
the same magnitude have the same statistical weight. On 
the other hand, by comparison with the fluctuation theo- 
rem (l4"9")l . it follows that here it is not possible to associate 
a temperature to the stochastic Levy dynamics. 

Other distinctive features of the problem can be char- 
acterized by analyzing the statistics of the dimensionless 
work (|3"4")1 . After a simple changes of variables, from Eq. 
(I5T)1) , the Fourier transform of its probability p T (w) reads 



(66a) 
(66b) 



G W (X) = (exp[iAuv]), 

= exp(i\wo — B T \X\ a ) 

where the coefficients are 

B T 



w a = I, 



B T = 



(67) 



From Eqs. (|33|) and ([BTjl we get the asymptotic behavior 
b fodt\v*(t)r 



B T 



r |Jo T <%*(*)]T' 



(68) 
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FIG. 2: (Color Online) (a) Levy probability distribution 
Pt(u>) denned by Fourier transform (|66[1 with a — 1.5. The 
inset details the peak region, (b) Plot of z q {ln 9 [p T (to)] — 
\n q \p T (—w)]} for each time r. The inset shows the collapse 
of these curves when plotting the normalized FR symmetry 
$0) = Z ;{\n q [p T (w)] - ln q \p T (-w)]}/r e , Eq. CD}. The po- 
tential velocity is constant, a c = 0, Eq. (|41fl . Hence, 5 = 2/3, 
Eq. I|7ip . In natural units (see text) the parameters are 
v* = 7 = k = 1. The noise intensity reads 6(711*) 2 = 2. 




FIG. 3: (Color Online) (a) Levy probability distribution 
p T {w) defined by the Fourier transform (166[) with a — 0.75. 
The inset details the peak region, (b) Plot of 2j"{ln g [p T (w)] — 
ln g [p r (— to)]} for each time r. The inset corresponds to 
$(tu) ee z^{ln q [p T (w)} - ln q [p T (-w)]}/T S , Eq. (70}. The 
potential velocity is constant, a c = 0, Eq. (1411) . Hence, 
S — —2/3 < 0, Eq. (|71[) . In natural units (see text) the 
parameters are v, =7 = ^ = 1. The noise intensity reads 
6(7t;») 2 = 2. 



For the velocity dependence (HIT) , it follows 



that regime, it satisfies the relation 



= (69) 



where d a = ^X+^'o' • Therefore, when a > 1/(1 + a c ) 
the width of the distribution diminishes, while for a < 
1/(1 + a c ) it increases with time. The former behav- 
ior is consistent with a LDP [see Eq. (|A1[) ]. Therefore, 
it should be possible to establishes a (generalized) FR 
symmetry for p T (w). In the second case, the typical fluc- 
tuations of the scaled work w T increase in size at higher 
times. This anomalous behavior is inconsistent with both 
a LDP and the law of large numbers [23] . Therefore, we 
expect that p T (w) does not fulfill any FR in this case. On 
the other hand, by analyzing p T (w) through a standard 
FR, the transition a J* 1 (a c — 0) leads to the different 
characteristic behaviors found in Ref. 1231 . 



For I to— wq\ > b\I g , the probability distribution p T (to) 
behaves as p T {w) ~ c a B T /\w — wq\ 1+<t [see Eq. (1521) ]. In 



-{ln q [p T (w)} - hig[p r (-to)]} oc a p hW, (70) 



where the coefficients read 



2(1 



ap h Kid' a b -(7W*) 2 , (71) 



with d' - J 2/a 



l rT (c' cr da' T ) 1 , and q is defined by Eq. IjM)). 
that is, q = ^if • On the other hand, the expression for 
a p h is only approximated because it is based on a map- 
ping with a q-Gaussian distribution. In general, the value 
of Zq , Eq. ([54)) , differs from that corresponding to a Levy 
distribution. 

Eq. (170)) is valid for any value of a £ (0, 2) if to per- 
tains to the power-law domain. Equivalently, it does not 
applies for [to — tool ^ B^ a , Therefore, when B T decrease 
(increase) in time the FR symmetry is valid (not valid) 
for almost any value of to. In fact, the transition in the 
behavior of the characteristic width B T , Eq. (|f)9")) . also 
determines when a LDP applies (5 > 0) or not (5 < 0). 
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From Eq. flTJ) it follows S ^ if a ^ 1/(1 + a c ). These 
properties are explicitly shown in the next figures. 

In Fig. 2(a) we plot the (exact numeric) Levy distribu- 
tion p T (w) |37| obtained from its Fourier transform (|66p . 
We assumed a constant potential's velocity, a c = in 
Eq. gl]), and a = 1.5. Hence, q = 1.8 [Eq. ([Hi])], and 
<5 = 2/3 > [Eq. (I7T1) ]. Consistently, the peaks around 
w = 1 diminish their wide for higher times. By using 
natural units, the unique free parameter is the (dimen- 
sionless) noise intensity [Eq. ([58])]. &To(7 2 v*/fc)~ cr . 

In Fig. 2(b) we plot the dependence with w of 
Zq{\n q [p T (w)} — ln q [p T (— w)}} for different times r. Small 
deviations with respect to a linear behavior are observed 
around |w;o — B^ IJ \ < w < \wq + B^' j \. Their magnitude 
diminish with time. In the inset we show the collapse to 
a single curve when introducing the time normalization 
factor (l/r 5 ), Eq. We checked that the same prop- 

erty is valid for higher times, indicating the consistence 
between the generalized FR and its associated LDP. The 
value of (XphTQ was estimated from the slope of the col- 
lapsed curves (inset). The theoretical estimation, Eq. 
Cfl]), gives a ph T$ = 1.85. 

In Fig. 3(a) we plot the Levy distribution p T (w) for 
a = 0.75, and a c = 0. Hence, q = 15/7 ~ 2.14 [Eq. (|64])]. 
and 5 = -2/3<0 [Eq. ([71])] ■ A negative S implies that 
the wide of the peaks around w = 1 grows with time. The 
linear behavior of z q {ln q [p T (w)} — \a. q [p T (— w)}} with w is 

only valid for w > w^ + B^ a , where B T increases in time. 
Added to this failure, after introducing the normalization 
factor (1/t s ), the curves does not collapse into a single 
curve (inset). These properties are parallel to the inap- 
plicability of a LDP. On the other hand, the theoretical 
estimation, Eq. (ITTI) . gives a p hTQ = 0.21. 

IV. SUMMARY AND CONCLUSIONS 



standard case when written in term of the LDFs, Eqs. 
(IA251) and (lA26l) . 

The general formalism was applied for characterizing 
the fluctuations of the work performed on a dragged par- 
ticle immersed in a complex environment. When the 
power-law nature of the dynamics is induced by (spa- 
ciotemporal) temperature fluctuations, the work statis- 
tics is given by a q-Gaussian distribution. The FR sym- 
metry is scaled by the environment average temperature. 
This novel fluctuation theorem [Eq. (|49[) ] may in princi- 
ple be confirmed in different experimental setups [27] . 

By analyzing the case in which the environment is rep- 
resented by an external Levy noise, we reinterpreted the 
results of Ref . [23| . Taking into account the superstatis- 
tical model, we conclude that some of those results are 
not valid in general. Due to the interplay between the 
noise statistics, the velocity of the power input, and the 
particle dissipative dynamics, in the long time regime the 
probabilities of positive and negative fluctuations of equal 
magnitude become identical. This result follows from the 
asymptotic vanishing of the characteristic constant that 
defines the work probability symmetry, Eqs. (|63[) and 
(1551 . On the other hand, the time-scaled work only sat- 
isfies the generalized FR when its behavior is compatible 
with a LDP, that is, the size of its characteristic fluctua- 
tions must to diminish with time (Fig. 2 and 3). 

While the uniqueness of the present proposal was 
not proved, based on the requirements that it satisfies 
[38}, one can conclude that it may be considered as a 
valid and solid tool for analyzing nonequilibrium fluctu- 
ations in systems characterized by power-law distribu- 
tions. The stretched relation with nonextensive thermo- 
dynamics [29| , as well as its applicability in specific exper- 
imental setups [36] are open problems that with certainty 
deserve extra analysis. 



We have introduced an alternative definition of FR 
symmetry that satisfies two conditions. In the regime 
where the probability of interest develops power-law tails 
the symmetry is expressed through a linear behavior. 
Furthermore, the generalized symmetry has associated 
a large deviation-like theory. 

The FR symmetry can be written as a difference be- 
tween the generalized q-logarithm of the probability dis- 
tributions for positive and negative fluctuations, Eq. ©. 
The parameter q depends on the exponent of the power- 
law tails. In terms of a generalized characteristic func- 
tion, Eq. (JSJ) , the proposed FR can be expressed as in the 
standard case, Eq. ([TU[). Similar relations, Eqs. ([2"Tjl and 
(1251) . were formulated for time-scaled variables. Based 
on large deviation-like principle, in the Appendixes we 
showed that a set of LDFs can be consistently defined 
for both the probability distribution and its associated 
characteristic function. The standard Legendre structure 
connecting them remains valid even in presence of self- 
similar power-law distributions, Eqs. (IA23[) and (| A24I) . 
Therefore, the generalized FR can be expressed as in the 
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Appendix A: Large deviation functions 

In this Section it is shown the consistence of the pro- 
posed FR with a large deviation-like principle. It is not 
obvious that an arbitrary generalized FR may satisfies 
this condition. Specifically, we show that it is possible 
to define two LDFs from the (long time) asymptotic be- 
havior of the probability density and its associated q- 
characteristic function. Both of them become related by 
a Legendre-Fenchel transformation [2(|. These results 
provide a solid mathematical support to the proposed 
FR. 
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Unnormalized scheme 



jointly with the inverse equation 



We base our analysis on time-scaled variables, Eq. 
(I20[) . A large deviation- like principle relies on providing 
a general structure for the probability p(/x) in the long 
time regime. Instead of a standard exponential structure 
poj , here we assume 

K^x^expJ-^-VM)]- (Al) 

As before, the symbol x denotes an equality valid in a 
long time regime. The factor t s (S > 0) measures the 
time-speed of ip(/j,) > 0. As in Refs. [2J, [25| (see also 
Appendix D of Ref. |20|), we consider the case in which 
6 =/= 1. On the other hand, the factor (t s ^ 2 /C q ) in front 
of the q-exponential is necessary for providing the rights 
units and normalization of p(fJ>). 

After a simple manipulation without involving any ex- 
tra approximation, Eq. (|Al|) can be rewritten as 



t 6/2 

p(aO ~—®q expJ-iVO)] 



where the exponent rj reads 



V 



(A2) 



(A3) 



Written is this way, given that lim t _ ! . 00 ln q [t s / 2 /Cql/t 11 
0, the function ip(n) can be obtained as 



-1 



ip(li) = Jim — \n q [p(p) 



(A4) 



Hence, it can be read as the probability's LDF [20(. 

Another LDF can be defined from the asymptotic time 
behavior of Z q {\). Its structure can be obtained from the 
definition (|2"3")l . after taking into account the probability 
asymptotic behavior (|A1[) . In general, the resulting in- 
tegral cannot be obtained exactly. Nevertheless, it can 
be worked out through a steepest descent approximation. 
In Appendix B we derive the asymptotic expression 



Z q (X) x exp= 



1 



where the index of the q-exponential reads 
- l + « 



<1 = 



3-q' 



(A5) 



(A6) 



Hence, the LDF associated to Z q (X) can be defined as 
-1 2 



cy\n q [z q (\)]. 



(A7) 



The steepest descent approximation establishes a link 
between both LDF (Appendix B). They are related by 
the Legendre-Fenchel transformation 



6(A) = min[(^(/x) + A/x], 



(A8) 



ip(n) — max[G(A) — A/x]. 



(A9) 



These relations also arise from a standard LDP, where the 
asymptotic behavior of the probability and its character- 
istic function scales with standard exponential functions. 
Remarkably, this Legendre structure remains valid even 
when the distributions develop power-law tails. 



a. Symmetries of the LDFs 

After establishing a large deviation-like principle [Eqs. 
Oj) and (| A5[) ] . we ask about the symmetries that the 
LDFs must to satisfy when the generalized FR is valid 
in the long time regime. A probability p(fi), with the 
asymptotic structure (|A2|) . satisfies the FR (f2"Tj) if the 
LDF (f(pJ) fulfill the condition 



(A10) 



Furthermore, the characteristic function (|A5|) satisfies 
the symmetry (|2"4"1) if the LDF 0(A) satisfies 



9(A) = 9(-A + a). 



(All) 



Both conditions are consistent between them. In fact, 
one can be derived from the other by using the Legendre 
structure defined by Eqs. (|A8[) and (| A9|) . The demon- 
stration is exactly the same than in the standard case 

ana. 



b. q- Gaussian distribution 

It is very instructive to exemplify the previous re- 
sults with an arbitrary q-Gaussian distributed variable. 
Let consider a stochastic variable x s t(t) whose long time 
statistics is given by Eq. (|13j) under the replacements 
xq — > Xo(t) and j3 — > f3(t). Furthermore, we assume that 
asymptotically these objects behave as 



x (t) x n t K , 



Ao 
t K ' ' 



(A12) 



Both, k and k' are positive exponents, while /io and Ao 
are characteristic constants. Notice that both the average 
(strictly the most probable value) and the characteristic 
width of the distribution p(x) grows with time. 

By defining fi st (t) = x st {t)/t K [Eq. J2DJ], using the 
change of measures p(^i)dfi = p(x)dx, from Eq. (| 13[) it 
follows the distribution 



^ y-jf ex P? [-A t 2K - K ( M - Mo) 2 ]- (A13) 

Hence, by comparing with Eq. (|A1|) it follow the identi- 
fications 



S = 2k- k' 



C„ 



N q 



(A14) 
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The probability LDF reads 



(Ao) ; 



1-q 



(A15) 



In order to be consistent with a LDP, the exponent 5 must 
be positive, 2k > k' . Hence, the width of the distribution 
p(n) diminishes with time. This property is expected for 
a time-scaled variable, Eq. (|20l) . A common situation 
corresponds to k = k' = 1 giving 5=1. Even in this 
case, the exponent rj, Eq. (|A3[) . which defines the limit 
(|A4|> . is different from one for q ^ 1. 

The q-characteristic function (1231) associated to the dis- 
tribution (|A13[) can be obtained exactly by using the pre- 
vious result (|17p. After a simple change of variables, we 
get 



Z 9 (A)xex Pg -{(^)t 5 C^ 



r X 2 AQ~ 9 



4(A ) 



A/x 



(A16) 



where g is given by Eq. (|T5J|, 5 and C g by Eq. (|A14|) . We 
notice that this structure corresponds to that obtained 
from a steeps descent integration, Eq. (|A5I) . In fact, that 
approximation is exact for a q-Gaussian distribution. 

By comparing Eq. (|A5[> and (| A16|) , we obtain the LDF 
9(A). It reads 

9(A) = - « _„ +A Mo . (A17) 
4(A ) — 

It is straightforward to prove that </?(//) [Eq. (|A15[) ] 
and 0(A) [Eq. (|A17[) ] are related by the Legendre- 
Fenchel transformations (|A8[) and (|A9j) . Furthermore, 
both LDF satisfy respectively the symmetries (|A10|) and 
(|ATT]) with the same constant a, which reads 



(Ao) ; 



(A18) 



As expected, when q 1 the expressions (|A15I) . (|A17[) . 
and (|A18|) reduce to those corresponding to a normal 
Gaussian distribution. 



Normalized Scheme 



Eqs. (|A10[) and (|A11|) are equivalent to the unnor- 
malized FR ([2"T]). The normalized FR ([2SJ can also be 
expressed through a set of renormalized LDFs. We define 

^(/i) - t lim |°° d/^)]*) Kb(M)], (A19) 

and similarly 



e""(A ph ) = t Um ^K[zf (X ph )\. 



(A20) 



In the asymptotic regime, the relation (1261) is equivalent 
to X ph = Alimt^oo J^oo dn\p(n)] q . In fact, it is pos- 
sible to demonstrate that 



lim 



1 



t s-r, 



q 



(A21) 

where (5 — rj) = (q — 1)5/2. Hence, rj is the same 
exponent defined in Eq. (| A3|) . The last equality in 
Eq. (|A21I) follows by calculating the integral through 
a steepest descent approximation, where p(u) is given by 
(|Aip . In the derivation we used the result (|B9|) and the 
e qua lity ^ = r(^)r(^)/[r(^)r(^)]. By 
comparing the LDFs corresponding to the unnormalized 
[Eqs. (TA4l and ([AT]) ] and normalized [ (|A19|) and (|A"20j) ] 
schemes, from Eq. (|A21|) it follows the relations 



@ ph {X 



ph) 



9(A P ^^- 1 



(A22) 



Taking into account Legendre-Fenchel transformations 
(|A8]) and fA9]), Eq. ([A22]) implies that 

e ph (X ph ) = min[^(/i) + Vm], (A23) 

jointly with the inverse relation 

if ph {n) = m^[e ph (X PJ ) - X phl j] . (A24) 



Therefore, the normalized definitions (| A19|) and (|A20j) 
also maintain the Legendre structure associated to a large 
deviation theory. 



a. Symmetries of the LDFs 



The (unnormalized) symmetries (|A10[) and (jAllI) . 
added to the relations (|A22j) . lead to the equivalent rela- 
tions 

- <fP h (n) + tp ph (-n) = a ph fi, (A25) 
and consistently 

Q ph (X ph ) = e ph (-X ph +a ph ), (A26) 
where the constant a p h reads 

rv 

(A27) 



Otph 



u 1 



b. q-Gaussian distribution 

For the previous example, defined by the q-Gaussian 
distribution (jAl3|t . it follows the LDF 



<^ fe (M)=fM) A o(M-Mo) 2 - 



(A28) 
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The q-characteristic function (|A16I) leads to 

A 2 

Zf(X ph ) x exp, {t s [ - \ phfM3 ] }. (A29) 



Therefore, we get 

o ph (x ph ) 



\2 



4(V) A o 



+ Xph^O- 



The constant 



reads 



ot p h = 4/i A ( 



3-q 



(A30) 



(A31) 



Notice that all dependences of the LDFs with the param- 
eter q can be absorbed in the coefficient Aq. 



Appendix B: Steepest descent approximation 

Here we develop a set of approximations that allow 
to calculate the long time behavior of Z q (X) from the 
asymptotic behavior of p(p). The procedure is similar to 
that given in Refs. [l(| (see Appendix C) and [20j | . 

By introducing (jAll) in Eq. (|2"3")l . after some calcula- 
tions steps, we get 

i<5/2 poo 

ZqW * yr / dnexp q {-lfCl-*[<p(n) + A M ]}. (Bl) 

This integral cannot be performed in an exact way. In 
order to proceed, we introduce an integral representation 
of the q-exponential function (q > 1) [23. |35| 



i 



r(^r) Jo 



dTT |5f e -r[l + (,-l)»] 



(B2) 



(2 > 0). Hence, after inverting the order of the integrals, 
we write 



Z g (X) 



t 5/2 



n^)c q Jo 



2-q _ 

d TT g-ie T 



dfx (B3) 



x exp{-t 5 r[^( M ) + A/*]Cj-«(g - 1)}. 



At long times t, the integral in the fi variable can be 
worked out with a steepest descent integration method. 
The main contribution to the integral comes from the 
value of /i, /i*, that minimizes the exponential. Defin- 
ing /j,* = /J* (A) by the condition <fi{n*) = — X, [<p(z) = 
(d/dz)ip(z)], we can approximate 

(pfa) + A M « (p(jf) + A M * + 2^(M*)(M ~ A**) 2 - (B4) 

By assuming that (p((J>) is a convex function to have 
a minimum, <p(fJ,*) > 0, and using the integral 



J^dfie^ 2 = (tt/o) 1 / 2 , it follows 



(B5) 



, dTT (^-h) e {-r[i+t s e(x)cl-"( q -i)]}^ 
'0 

Here, 8(A) denote the function 

9(A) = <p(n) + A M , (B6) 
where the value of fi follows from the condition 

0(M) = -A. (B7) 



By using again the integral representation (|B2I) , we arrive 
to the expression 



Z q (X) x A q exp= 



3-9 



(B8) 



constant is 



where q is given by Eq. (|A6|) . g = (1 + g)/(3 — g). The 



4, = - 



r(^ T )c 9 2 Vfe-i)^*) c q - 

(B9) 

Therefore, Eq. (|B8[) leads to the expression (|A5I) . In 
the previous expression, the second equality follows from 
Eq. ([T4|) . while the last estimation follows by apply- 
ing a steepest descent approximation to the condition 
f^°° dfip(u) x 1, where p(u) is given by (|Alj) . jointly 
with the condition that the LDF ip(fi) vanishes at its 
minimum Mmin, ^(Mmin) = 0, which consistently implies 
9(0) = 0, Z q (0) x 1. 

Eqs. (|B6| and ([B7]) show that 9(A) is the Legendre 
transform of <p(/i), which in turn can be written as in Eq. 
(|A8|) . 9n the other hand, the derivative of 9(A) with 
respect to A is 



which from Eq. (|B7p leads to 

9(A) = M . 



(B10) 



(Bll) 



This shows that 93 (m) is given by the inverse Legendre 
transform of 9(A), 



ip{y) = 9(A) - Am- 



(B12) 



Here the value of A follows from the condition (jBllI) . 
leading to Eq. (IA9|) . In fact, by taking the derivative of 
(IB11[) with respect to A and using the derivative of (|B7|) 
with respect to m we can confirm that ip(fx) is convex 
because 9(A) is concave, </3(m) = — 1/6(A). 
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